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We present a theoretical model and respective numerical simulation of 
cascaded four-wave mixing of femtosecond pulses in bulk x media, ev- 
idencing the importance of two-dimensional interaction geometries in the 
efficient generation of frequency-converted beams, and in agreement with the 
experimental results in terms spectral and spatial distribution and efficiency, 
demonstrating compression without the need of complex amplitude or phase 
control. © 2009 Optical Society of America 

OCIS codes: 190.4380, 190.7110, 190.3270, 190.2055. 

The simultaneous phase-matching of multiple nonlinear optical processes can give rise to 
interesting cascading phenomena, with each consecutive process efficiently making use of op- 
tical signals generated by the previous ones. In particular, cascaded second-order processes [1] 
based on sum- and difference-frequency generation of intense laser pulses in x media can 
be very useful for generating new frequencies [2] beyond the usual signal and idler waves 
produced in single-step processes. However, these can only occur in non-centrosymmetric 
materials. Third-order nonlinear processes, namely four- wave mixing (FWM), have the ad- 
vantage that x is nonzero for all media, without any restrictions imposed by symmetry. 
Therefore, FWM can occur in any material, including isotropic media such as glasses, gases 
and plasmas. In particular, nondegenerate FWM of noncollinear laser pulses with frequencies 
Uq and uoi (loi > u ) can be geometrically phase- matched in media with positive dispersion 
to generate a third order difference-frequency signal u>2 = 2uj\ — ujq [3]. Cascading of this 



process (generating pulses at u n = nui — {n — 1)ujq where n is the beam number or order, 
with negative values denoting frequency-downshifted beams) would seem at first impossible 
due to the need of different phase-matching angles for each intermediate step. However, this 
condition can be greatly relaxed by using broad bandwidth pump pulses. 

In the first demonstration of highly nondegenerate cascaded four-wave mixing (CFWM) 
of femtosecond pulses in bulk x media [4] it was shown that two noncollinear pump pulses 
with different frequencies overlapping in a nonresonant Kerr medium can generate a series of 
angularly separated multicolored pulses with broad bandwidths and frequencies extending 
from the infrared to the ultraviolet. Broadband phase- matching of the corresponding multi- 
ple processes was achieved by adjusting the interaction angle between the pumps, and the 
use of a thin slide as the nonlinear medium minimized material dispersion and competing 
nonlinear optical effects such as self- and cross-phase modulation while maintaining high 
conversion efficiencies (approx. 10%). The efficient extension of CFWM generation down 
to 160 nm in gas-filled hollow waveguides [5], with phase-matching achieved by tuning the 
gas pressure, further evidenced the universal nature of nonresonant CFWM processes. More 
recently, CFWM was also observed in noncollinear four-wave parametric amplification of 
ultrashort laser pulses in water [6] . Another interesting characteristic of nonresonant FWM 
and CFWM processes is the possibility of generating and synthesizing compressed pulses. 
This is again a direct consequence of the modulation of the incident fields by the material 
excitation, in its turn proportional to the product of the initial pump field envelopes. Fuji et 
al. obtained ultraviolet 266 nm pulses with 12 fs in duration by FWM of the fundamental 
and second harmonic pulses from a 25 fs Thsapphire laser in a gas filament, where cascaded 
generation of radiation at 200 nm was also observed [7]. In recent work, polarization gated 
cross-correlation frequency-resolved optical gating (PG-XFROG) measurements of several 
simultaneous cascaded beams generated in fused silica and covering 1.5-octaves in band- 
width were performed that support the possibility of using CFWM to Fourier-synthesize 
high-power optical pulses in the single-cycle regime [8,9]. 

The interest and characteristics of CFWM processes have been calling for theoretical and 
numerical models capable of accurately describing its detailed features and potential. Hart 
et al. [10] used coupled wave equations as well as the one dimensional nonlinear Schrodinger 
equation to describe multiple FWM of nanosecond pulses in optical fibers. Lichtman et 
al. [11] obtained exact analytical solutions of the total field for ID propagation of 2 or 
3 monochromatic input waves but assuming a nondispersive medium and perfect phase 
matching. 

We found out that in order to accurately reproduce the main features of broadband CFWM 
and to account for geometrical phase matching, a two-dimensional approach is required. In 
this letter we present a numerical study of CFWM of noncollinear femtosecond pulses with 
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Gaussian spatial profiles crossing in a thin fused silica slide. In our model we consider the 
propagation of the pulse envelope on the 2D plane defined by the pump beams, using a 
reference frame that moves along the propagation axis z with the group velocity of one of 
the pumps, as given by Eqs. 1-2. 
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where A is the envelope of the pulse with electric field E (x, t) = A (x, t) exp (ik z — iuj t), 
5uj = uj—ujq, v g is the group velocity at the reference frequency u , k (u) = n^ui/c, k = k (u ) 
and k' (u) = ko + 5uo/v g . Eq. 1 describes linear propagation in the frequency domain, and 
includes diffraction and dispersion (although due to the small longitudinal dimension of the 
medium the group velocity mismatch and pump walk-off are practically negligible). Eq. 2 
contains the nonlinear terms in the time domain, which comprise the instantaneous nonres- 
onant Kerr effect and self-steepening. The accurate evaluation of the self-steepening terms 
is essential because the generated spectral width is of the same order as u Q . No additional 
delayed response terms were required to faithfully reproduce the experimental results. 

We solve these equations in a split step fashion, the first using a (2,2)-Pade approximant for 
wide-angle propagation as two (1,1) steps [12], and the second using a second order Runge- 
Kutta method, with time derivatives approximated by 5-point finite differences. We used 
the appropriate parameters to simulate our experiment described in [9]. The green (565 nm, 
60 fs, 1.5 x 10 12 W/cm 2 ) and orange (614.8 nm, 80 fs, 1.8 x 10 12 W/cm 2 ) transformed-limited 
pump pulses propagate in 172 /jm of fused silica (corresponding to a 150 /xm thick slide tilted 
at 45 degrees, which reduces the angle between the normal to the slide and the generated 
frequency upshifted beams thus preventing total internal reflection at the silica-air interface), 
crossing at an internal angle of 1.57 degrees. The green pulse is also delayed 48 fs with respect 
to the orange at the entrance of the slide, which resulted in better defined spectra for the 
CFWM orders, in agreement with the experimental measurements. 

The simulations provide the time dependent field distribution at the end of the slide as 
a function of the transverse coordinate x and frequency uj. The corresponding spectrum 
(obtained by integration along the transverse direction) is given in Fig. 1, clearly showing 
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the asymmetric generation of up and down-converted orders. The efficiency of the generated 
frequency upshifted pulses follows the same pattern as the experiment, although it has a 
sharper cutoff after the 12th order that depends on the interaction angle and the pulse 
energy This discrepancy is possibly due to experimental calibration errors which are larger 
at the spectral edges. For collinear pumps the efficiency is low because the process is not phase 
matched. As the angle is increased (> 1.4 degrees) phase matching favors the generation of 
higher order modes, as shown in [9]. At larger internal angles (> 2 degrees) the overall 
efficiency drops rapidly. 

The 9 — A spectrum (Fig. 2) shows the main features of the CFWM process in the cor- 
responding domains. The various orders are emitted with different wavelengths and angles, 
with decreasing angular separation between beams of increasing order. 

The simultaneous spatial and temporally superimposed CFWM orders have been charac- 
terized experimentally with PG-XFROG where a portion of the orange pump pulse acted 
gate [8,9]. Here, to simulate the focusing and obtain the total electric field needed to calculate 
the PG-XFROG traces we simply subtracted the linear phase (group delay) acquired due to 
the transverse propagation through the fused silica slide. The resulting field intensity and 
respective spectrum are shown on Fig. 3 which shows a train compressed pulses separated by 
the pump beat period (23 fs) as expected. The central peak has a duration of 3.1 fs FWHM. 

Fig. 4 compares the experimental and simulated PG-XFROG traces for the second 
frequency-upshifted pulse and the first downshifted pulse. We see that the simulation results 
are in good agreement with the measurements, providing similar values for the bandwidth 
and duration of the generated pulses (which are also shorter than the pumps). The observed 
difference in the central frequency of the upshifted pulse is due to the slight changes in the 
central frequency of the pumps that can occur between experiments. The pulse also shows 
some positive chirp due to the nonlinear phase imposed by the pumps. This is more clear 
in Fig. 5, which shows the calculated PG-XFROG traces for the synthesized field of Fig. 3 
assuming gate pulses of 20 and 80 fs, that provide more detailed information in the temporal 
and spectral domains, respectively [8,9]. The right trace (80 fs gate) reproduces the condi- 
tions of [8,9], and is a good match to the experiment. The shorter gate pulse should enable the 
direct observation of fine temporal structure of the synthesized pulse train, provided that the 
relative phase between the pump pulses is kept constant throughout the measurements [9]. 

In conclusion, we developed a theoretical model capable of accurately describing the main 
features of highly nondegenerate CFWM of noncollinear femtosecond pulses in the spatial, 
spectral and temporal domains. A nonlinear propagation equation in two spatial dimensions 
was introduced, without imposing artificial restrictions on either the number or frequencies of 
the generated orders, and presented a numerical method to solve it efficiently. The simulations 
show that the instantaneous Kerr response and self-steepening are necessary and sufficient 
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to reproduce the experiments and demonstrate compression without complex amplitude 
or phase control. They may provide a tool to study new configurations and materials for 
generating and compressing pulses in previously hard to access domains. 

J.L.S. acknowledges financial support under Grant No. SFRH/BD/24523/2005 from FCT- 
MCTES, Portugal. 
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List of Figures 

1 Spectrum of the multiple CFWM orders generated in fused silica. The spec- 
trum spans over 1.5 octaves (the two pump frequencies correspond to the two 
larger peaks). Dashed line: Experimental results (see [9]), Solid line: Numerical 
simulation, Dotted line: Numerical simulation for collinear pumps 8 

2 (Color online) Simulated 9 — A spectrum of the CFWM beams, clipped at 10~ 6 

of the maximum 9 

3 Recombined and focused CFWM pulses, (a) Electric field as a function of the 
delay (in fs). (b) Intensity, (c) Spectrum (logarithmic scale) 10 

4 (Color online) Experimental (measured and retrieved) and simulated PG- 
XFROG traces of the second frequency-upshifted pulse (top row) and first 
downshifted pulse (bottom row). The gate pulse is the 80 fs input orange 
pump. The retrieved pulse durations are 30.7 and 40.2 fs for second upshift and 
first downshift respectively (FROG error was 0.014 and 0.008 for a 128x128 
grid) 11 

5 (Color online) Simulated PG-XFROG, clipped at 10 -4 of the maximum, for 
the recombined and focused CFWM field assuming a 20 fs gate pulse. Right: 

80 fs gate pulse 12 
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Fig. 1. Spectrum of the multiple CFWM orders generated in fused silica. The spectrum spans 
over 1.5 octaves (the two pump frequencies correspond to the two larger peaks). Dashed 
line: Experimental results (see [9]), Solid line: Numerical simulation, Dotted line: Numerical 
simulation for collinear pumps. 



8 



"3) 
d 

3 
"3 

d 



-10 




I 



400 500 600 700 800 
Wavelength (nm) 



r° 

-6 
-12 

-18 
-24 

-30 « 

-36 
-42 
-48 
-54 
-GO 



Fig. 2. (Color online) Simulated 9 — A spectrum of the CFWM beams, clipped at 10 6 of the 
maximum. 
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Fig. 3. Recombined and focused CFWM pulses, (a) Electric field as a function of the delay 
(in fs). (b) Intensity, (c) Spectrum (logarithmic scale). 
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Fig. 4. (Color online) Experimental (measured and retrieved) and simulated PG-XFROG 
traces of the second frequency-upshifted pulse (top row) and first downshifted pulse (bottom 
row). The gate pulse is the 80 fs input orange pump. The retrieved pulse durations are 30.7 
and 40.2 fs for second upshift and first downshift respectively (FROG error was 0.014 and 
0.008 for a 128x128 grid). 
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recombined and focused CFWM field assuming a 20 fs gate pulse. Right: 80 fs gate pul: 
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